Full control by locally induced relaxation 
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We demonstrate a scheme for controlling a large quantum system by acting on a small subsystem 
only. The local control is mediated to the larger system by some fixed coupling Hamiltonian. The 
scheme allows to transfer arbitrary and unknown quantum states from a memory on the large system 
("upload access") as well as the inverse ("download access"). We study sufficient conditions of the 
coupling Hamiltonian and give lower bounds on the fidelities for downloading and uploading. 
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The unitarity of Quantum Mechanics implies that the 
information lost by an open system during its dynami- 
cal evolution must be contained in the environment and, 
possibly, in the correlations between the system and the 
environment [l| . In the context of repetitive applications 
of the same quantum transformation this fact has been 
exploited to achieve noise protection [2j], cooling, state 
preparation 0, Q| , and quantum state transfer [H] . The 
mathematical aspects related with the convergence of the 
associated trajectories are well studied (see Ref. @ and 
references therein). It has however been largely over- 
looked that the relaxation effects associated with these 
approaches can also imply full control by acting on a local 
subsystem only. This is analogous to what happens in the 
case of the universal quantum interface of Ref. [7[ where 
the local control is mediated to the whole system by some 
fixed coupling Hamiltonian. Once this is achieved, apart 
from cooling and state preparation, it is also possible to 
perform arbitrary quantum data processing (e.g. mea- 
surements, unitary rotations). 

In this paper we discuss an explicit protocol for uni- 
versal control of a composite system by operating on it 
with a simple repetitive local quantum transformation 
and we provide lower bounds for fidelities obtainable af- 
ter finitely many steps. Furthermore turning the prob- 
lem into a graph theoretic one we provide an easy-to- 
check sufficient criterion to verify if a given global net- 
work Hamiltonian is capable of mediating control. 

The results presented here pave the way to new appli- 
cations of quantum control and quantum computation. 
Arguably our downloading and uploading protocols (see 
below) may be useful for the control of quantum hard 
drives and quantum RAM Q, or CCD- like application 
for the external control of permanently coupled arrays of 
sensors (see caption of Fig [IJ . A large scale experimen- 
tal realization of the scheme discussed in this paper is 
not realistic at the moment since it requires the ability of 
performing full quantum computation on a large memory 
system. However proof-of-principle tests could be prob- 
ably realized in hybrid quantum networks by exploiting 
the methods proposed in Ref. (e.g. controlling a per- 
manently coupled optical lattice - where local control is 
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Figure 1: The systems C and C are coupled through the time 
independent Hamiltonian H . The system C is controlled by 
performing regular swap operations Si between it and a quan- 
tum memory Mi. In CCD-like application C is an array of 
permanently coupled sensors that are used to probe external 
signals (schematically sketched by the wiggling curve) : infor- 
mation from the sensor is extracted through C. 



generally difficult - by coupling it to a fully controllable 
array of trapped ions). 

Protocol:- We consider a tripartite finite dimensional 
Hilbert space 7i = He <8> Hq (g> %. Full control (i.e. the 
ability to prepare states and apply unitary transforma- 
tions) is assumed on system C and M, but no (direct) 
control is given on system C. System C and C are cou- 
pled by some time-independent Hamiltonian H. We will 
show that under certain assumptions, if the system CC is 
initialized in an arbitrary state we can transfer ( "down- 
load") this state into the system M by applying some 
simple operations which act locally on M and C. Like- 
wise, by initializing the system M in the correct state, we 
can "upload" arbitrary states on the system CC. These 
two schemes ensure full controllability of CC: for in- 
stance one can perform arbitrary quantum operation on 
such system by transferring its state into M, applying 
the equivalent operation there, and transferring the re- 
sulting state back to CC. In this context M functions as 
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a quantum memory and must be at least as large as the 
system CC. As sketched in Fig. [T] we can imagine it to 
be split into sectors Me, 7i M = Hm, with the Hm £ 

being isomorphic to He (i-e. dimHM e = diraHc)- 

For downloading we assume that the memory M is 
initialized in \e)M = 0^ \ e )M e with the vectors |e) to be 
defined in the following. To download an arbitrary initial 
state \tp)cc of CC into the memory M we perform a 
sequence of unitary gates between M and C, intermitted 
by the time evolution U — exp[— iHt] on CC for some 
fixed time interval t. More specifically, at step i of the 
protocol we perform a unitary swap St between system C 
and system Mg. The protocol stops after the Lth swap 
operation. The resulting global transformation is thus 
represented by the unitary operator 



W = USlUSl^x ■■■USv US!. 



(1) 



As we will see in the next section, the reduced evolution 
of the system C under the protocol can be expressed in 
terms of the completely positive trace preserving (CPT) 
map r defined in Eq. ([4]). Our main assumption is that 
the system C is relaxing [l2| under repetitive application 
of t, i.e. linin^oo r n {p) — p* for all initial states p. This 
behavior is also called mixing 2, ll[ or absorbing Q . In 
what follows we will focus on the case in which is a 
pure state \E)q{E\. When this happens it is possible to 
show that, for sufficiently large L, the transfer of \ip) from 
CC into M can be done with arbitrarily high fidelity [|| 
and the transformation which allows one to recover 
from M can be explicitly constructed. 

For uploading an arbitrary input state from M to 
CC one is tempted to revert the downloading protocol. 
Roughly speaking, the idea is to initialize the memory in 
the state that it would have ended up in after applying 
W if system CC had started in the state we want to 
initialize. Then we apply the inverse of W given by 



W* = • ■ • S e U* ■ ■ ■ S L -iU ] S L tf 



(2) 



We will see that indeed this induces a unitary coding on 
M such that arbitrary and unknown states can be initial- 
ized on CC. The reader has probably noticed however 
that the transformation (|2|) is generally unphysical in the 
sense that it requires backward time evolution of CC, i.e. 
one has to wait negative time steps between the swaps. 
For this reason, even though the transformation origi- 
nated from is coherently defined at a mathematical 
level, it cannot be considered as a proper uploading al- 
gorithm for transferring states from M to CC: to stress 
this we will call the transformation associated to Eq. @ 
the reverse-downloading protocol. A proper uploading 
algorithm will be defined in the final part of the paper 
by imposing an extra hypothesis on the CC couplings 
and by adopting a simple change of perspective. For 
the moment we neglect this point and simply focus on 
the convergence properties of the downloading and the 



reverse-downloading algorithms associated with Eqs. (p} 
and H}. 

Cooling: - We start by showing that the action of 
W on CC is effectively equivalent to a cooling pro- 
cess which transfers any initial state into je^l-E^c-Let 
\ip)cG e ^cc be an arbitrary state. We notice that the 
C component of W\ij)) C c\^}M is always \e)c- Therefore 
we can write 

W\^) C c\e) M = \e) c [/v\E)c\4>)m + ^~v\^)cm] (3) 

with \A)q M being a normalized vector of CM which sat- 
isfies the identity q(E\A)q AI = 0. It is worth stressing 
that the decomposition ^ is unique and that rj, \4>}m 
and \A)cm are typically complicated functions of the in- 
put state |V , )cc- The quantity rj plays an important role: 
it gives us the fidelity between the initial state of CC and 
the target state \e)c\E)Q of the cooling process. An ex- 
pression for -q can be obtained by focusing on the reduced 
density matrix of the subsystem C. From our definitions 
it follows that after the first step of the protocol this is 

t( Pc ) = tr CM [USi (MccM ® |e) Af <e|)5 1 l7t] 

= tv c [U{pc®\e)c{e\)U^] , (4) 

with Pq = trc [|'0)cc(V'|] being the reduced density ma- 
trix associated with the initial state \ip)cc- Reiterating 
this expression we notice that the state of C after L steps 
can be obtained by successive application of the map ([!]). 
Consequently Eq. ([3]) gives r\ = ^{E\t l (pq) \E)q, which, 
according to the mixing properties of r given at the be- 
ginning of the section, shows that 77 — ^ 1 for L — > 00. 
Specifically we can use [l2| to claim that for all input 
states the following inequality holds 



\n-l\^\\r L (pc)-\E)c(E\h<KK L L d 



(5) 



where K is a constant which depends upon d(j = dimTYg 
and k G]0, 1[ is the second largest of the moduli of eigen- 
values of the map r. 

Coding transformation:- Let us now derive the de- 
coding/encoding transformation that relates states on 
the memory M to the states of CC. The idea is to ap- 
ply the decomposition |3]) to each element of a given or- 
thonormal basis {|^fc)cc} °f ^cCi ano - to define the lin- 
ear operator D on TLm which, for all k, performs the 
transformation 



D\tp 



k/M 



M 



(0) 



In this expression \ipk)M are orthonormal vectors of M 
used to represent the states \tpk)cc °f ^cc on M (for- 
mally they are obtained by a partial isometry from CC 
to M). The vectors \<pk)M instead are connected to the 
\ipk)cG through Eq. ([3|). Typically, for finite values of L, 
the \4>k)M will not be orthogonal. However it is possible 
to show that they become asymptotically orthogonal in 
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the limit of L 1. To see this we use the unitarity of 
the transformation W and the orthogonality of \ipk)cc- 
Indeed from Eq. ([3]) one can easily verify the following 
identity 



Skk' = y/rnVk' M{4>k\<f>k')M 



(J) 



+ -%)(!-%') CM(A fc |A fc ,)cM 



Defining rj = minfcryj,, we notice that for sufficiently 
large L this is a strictly positive quantity and converges to 
1 — see Eq. From the identity ([7]) it follows then that 
for k ^ k' one can write \M{<t>k\<t>k') m\ < (1— »7o)/?7o- This 
can now be used to bound the eigenvalues Xj of the linear 
operator D'D. Indeed the Cauchy-Schwartz inequality 
yields |Aj - 1| < d c0 (1 - rjo)/i]o, with d cC = dim7Y ce . 
Take now a polar decomposition D = PV with P being 
positive semidefinite. According to [lj, p 432] V is the 
best unitary approximation of D. In our case it satisfies 
the relations 



\D-V\\ 



i 



«S d cC (d cC -1)(1- Vo)/Vo 



(8) 



This is a key equation: thanks to Eq. (O, it shows that 
D can be approximated arbitrary well by the unitary 
operator V for L — > oo. 

Fidelities:- In what follows we will use V+ and V 
as our downloading and reverse-downloading transfor- 
mation, respectively. In particular, will be used to 
recover the input state |V>)cc °f the chain after we have 
(partially) transferred it into M through the unitary W 
(i.e. we first act on \i>)cc\ e )M with W, and then we 
apply V' on M). Vice- versa, in order to upload a state 
\"4>)m on CC by using the reverse-downloading protocol 
we first prepare CC in \eE) C Q then we apply to it 
the unitary transformation V and finally we apply W* . 

With these choices the fidelity for downloading a state 
\ip) M is simply given by F d (ip) = m(iP\V^ pm V\ip) M 
where pm is the state of the memory after W, i.e. 

p M ee tr ce [W(\i>) C cm ® |e)M<e|)W t ] 

= V \4>)m(4\ + (i-v) o-m , (9) 

(here we used Eq. ([3]) and introduced the density ma- 
trix om = trp[|A)(j M (A|]). Analogously the fidelity 
for transferring a state from M to CC through the 
reverse-downloading protocol is given by 

F rd (4>) ee ce <#r M [wV(|V>M<Vl 

A bound for F r and F rc ; follows by noticing that both 
these quantities satisfy the inequality 



For Fd this simply comes by expressing it in terms 
of Eq. ([9]) and by neglecting a positive contribution 
proportional to 1 — rj. For F rc i instead the inequal- 
ity (flT))) follows by replacing the trace over M with 
the expectation value on |0)m an d by using Eq. ([3]). 
One can now estimate the scalar product on the right 
hand side of Eq. (fTQ)) by observing that |at(0|V|^)m| ^ 
\ M {<t>\D\i>) M \ - \m(4\D - V\il>) M \. This can be further 
bounded by employing the inequality (jSJ) and the fact 
that \ M (<f>\D-V\i/j) M | < ||D-V|| 2 . If W)m isavectorof 
the basis \tpk)M, then \m(4>\D\i/j) m\ — 1 by the definition 
- see Eq. ([6]). For generic \i/>)m instead some simple alge- 
bra yields ^/rj \ M {<j)\D\^) M \ ^ y/rfo - d cC (1 - Vo)/Vo- 
Replacing all this into Eq. (fT0|) we finally get 



F > r]o - 4 d cC Vo)/Vo 



(11) 



which holds for F = F r) F r d- This is a lower bound for 
the fidelity of the downloading and reverse-downloading 
protocols: it is probably not tight but it is sufficient to 
show that F r and F r d converge to 1 in the limit of large 
L . According to Eq. ^ such convergence is ex- 
ponentially fast in L even though, not surprisingly, the 
bound deteriorates as the size d^ of the controlled system 
increases. 

Uploading protocol:- Let us now come back to the 
question about the operation W* being unphysical. To 
define a proper uploading protocol consider a modi- 
fied scenario in which CC is replaced by a isomorphic 
system CC 1 characterized by the Hamiltonian H' 

1 1 In this scenario the downloading protocol is de- 
scribed by the operator W = U 1 S L U' S L ^- ■ -U' S x 
with U' ee exp[— iH't] while the corresponding reverse- 
downloading transformation by the operator (W)t = 
Si(liy ■■■S L - 1 (U'^S L {U')1. Since U 1 = W it is 
not difficult to observe that the "unphysical" reverse- 
downloading algorithm of C'C induces a proper up- 
loading transformation for CC. To ensure that such 
algorithm converges it is hence sufficient to study the 
downloading and reverse-downloading protocols associ- 
ated with W and W. According to our previous results 
this can be done by focusing on the CPT map 



r'(p c )=tv c [17+ [p G ®\e)c{e\)U] 



(12) 



F(ip) > v\M(<t>\V\iP) M \ 



(10) 



which replaces r of Eq. It is then sufficient to as- 
sume t' to be ergodic with pure fixed point \E)q. When 
this happens we can define a transformation D' as in ([6|) 
and its unitary part V . The latter is the coding trans- 
formation which will be used for the uploading protocol 
of CC. Consequently the fidelity F up (ip) associated with 
such algorithm is bounded as in Eq. (fTTj) with rjo being 
lower bounded by Eq. ([5]) where the parameters K and k 
of t have being replaced by the corresponding quantities 
of t'. It should be noted that the definition of uploading 
protocol given here is more general than in other schemes 
relying on time-reversal symmetries 0] ■ 
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Figure 2: Examples of graphs coupled by Heisenberg-like in- 
teraction that can be controlled by acting on the black qubits 
only. The lower right graph provides a counterexample. 



A condition for controllability:- The mixing proper- 
ties of r and t' are typically independent (see for instance 
Ref. i): this makes it difficult to give a general condi- 
tion for the full controllability of CC (i.e. convergence of 
both the downloading and uploading protocols). Notably 
however a generic statement can be made using a result 
of Ref. " 



For the sake of simplicity here we will focus 



on the case in which CC is a network of coupled spins 
1/2 particles. According to II] we have that i) if the 
Hamiltonian H of CC preserves the number of spin exci- 
tations and ii) the vector |e)c|£')c7 is the only eigenstate 
with C in \e)c then the map r is mixing with fix point 
\E) C (here \e)c and \E)q represent states with all spins 
aligned down). However the Hamiltonian H 1 = —H as- 
sociated with t' has the same eigenvectors of H: hence 
the conditions i) and ii) also yields a sufficient criterion 
for determining that f is mixing with fixed point \E)(j. 

The analysis further simplifies by focusing on a two- 
sites interaction Hamiltonian. In this case a simple re- 
cursive analysis is sufficient to check if H satisfies the 
condition ii). For linear chains of spins this was dis- 
cussed in Ref. @: here we generalize this argument to 
arbitrary topology. To do so, define the graph Q with 
the spins of the network as vertices, and the non-Ising 
components of H as edges. Introduce also the following 
color code: a black vertex corresponds to a spin in | |) 
while a white vertex corresponds to a generic spin config- 
uration (i.e. not necessarily 1 |)). Consider now the case 
in which the subset C of the graph has all black vertexes. 
Our goal is to determine if such configuration is compat- 
ible with being a non trivial eigenstate of the network 
Hamiltonian (that is an eigenstate in which not all the 
vertices are black): if not, then the whole network can 
be controlled. This can be checked by noticing that the 
excitations (i.e. the white vertices) tend to propagate 
along the edges when H is applied to the graph: con- 
sequently only certain distributions of black and white 
vertices are compatible with the eigenvector structure of 
H (they must allow certain interference effects that pre- 
vents the propagations of the white vertices) . To exploit 
this property we introduce the following cellular graph 
automata: iff a vertex is black and has exactly one white 
neighbor, then this neighbor will turn black in the next 
step. Otherwise, vertices remain unchanged. It follows 
that the final state associated with a generic initial con- 



figuration is the one which is compatible with being an 
eigenstate of H and has the minimum allowed number of 
black vertices. Therefore to verify if the whole graph is 
controllable by operating on C it is sufficient to initialize 
C in the all-black configuration and C in the all-white 
configuration and let the graph evolve. If the final result 
has only black vertices then we conclude that C induces 
full control on CC (see Fig. [2] for examples of controllable 
graphs). Note that this method allows us to say some- 
thing about the structure and the asymptotic dynamics 
of Hamiltonians that are in general very far from analytic 
solvability, including disordered systems. 

Conclusion:- We have shown that a for large class of 
physically realistic Hamiltonians, a Quantum Computer 
can fully control a large system by inducing a relaxation 
on a small subsystem only. The fidelity converges expo- 
nentially fast with the size of the memory, so the required 
overhead of resources is low. The results discussed here 
can be generalized to maps r and r' are mixing with 
non-pure fixed point p» . In this case the fidelities of the 
protocols will be not necessarily optimal. Still one can 
provide non trivial lower bound for this quantities which 
depends upon the purity of p*. DB acknowledges the 
Swiss National Science Foundation (SNSF) for financial 
support. 
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